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SUMMARY 


A theoretical analysis of an airplane automatic control system 
incorporating a con 5 )ensating network is presented. The compensating 
network is a computing network that has characteristics that are the 
inverse of the airframe; consequently, airplane dynamics are eliminated 
from the system response. The transfer function of this type of control 
is developed, and the result is applied to the analysis of roll control 
systems. The basic roll control system is a displacement or bank-angle 
control system. This basic system was modified by feeding back rolling 
velocity and acceleration to the input of the compensating network to 
obtain velocity and acceleration control systems. In addition to the 
linear analysis, an analog-computer study was made to determine the 
effects of limiting the control-surface rate and displacement and of 
incomplete compensation on the command response of the roll control 
systems . 

The results are presented in the fonn of time histories of the 
lateral airplane variables and control-surface motions . 

For the forward-loop compensating network, the effects produced by 
limiting and incomplete compensation on the response characteristics of 
the system indicate that the system cannot at present be considered as 
a satisfactory automatic control system for interceptor airplanes. 


INTRODUCTION 


The manned all-weather interceptor has assumed an important role 
as an air defense weapon. Since it is anticipated that the interceptor 
may be automatically controlled in the attack phase of the intercept 
mission, an effective automatic control system is a necessary conponent 
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of the automatic interceptor system. The interceptor airframe, as 
presently conceived, is a monowing design that rolls to turn toward the 
interception point. Therefore, an automatic roll control system is a 
necessary part of the guidance system for these interceptors. As part 
of a general investigation of the attack phase of the automatic inter- 
ceptor, a roll control system incorporating a compensating network has 
been studied. 

Reference 1 proposes the use of a compensating network, as part of 
an automatic pilot, to compensate all or part of the lateral and longi- 
tudinal modes of the airframe. Compensating networks that wholly or 
partially cancel the longitudinal modes of the airframe, when controlled 
by a hioman pilot, have been discussed in references 2 and 5- 

The compensating network, as used in the roll control system of the 
present investigation, is a control-deflection computer designed to 
eliminate the lateral modes of motion from the response of the system 
to a command input. If the roll control system operates in a linear 
manner and if the con^jensating network is correctly designed, cancellation 
of the airframe dynamics is always obtained. However, in this investiga- 
tion physical limits were imposed on control- surface rate and displace- 
ment; thus the system became nonlinear. Also, small inaccuracies in the 
mass and aerodynamic data required to design the compensating network 
and changes in the airplane flight condition were used to introduce 
incomplete compensation into the system operation. The effect of these 
nonlinearities and inaccuracies on the behavior of the system in response 
to command inputs for a forward-loop compensating network is discussed in 
this paper . 

The results are presented as time histories of the lateral motions 
of the interceptor and control-s.urface motions, which were obtained on 
the Reeves Electronic Analog Computer (REAC) at Project Cyclone euad at 
the Langley Aeronautical laboratory. Some of these results have been 
previously summarized in reference 4 . 

r 

SYMB0I5 


t time, sec 

D differential operator, ^ 

a arbitrary constant in exponential input J^q(i - e”®"^) 

F(d) corapensating-network transfer function, forward loop 
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Fi(5) 

S(D) 

Si(D) 

G(D) 

N(D) 

6(D) 

Gi(D) 

n 

n(D), x(D) 
M 

0(D) 

T 

K 

R(D) 


Kl 

K2 

b 

S 

W 

g 

Kx 


compensating-network transfer fxmction, feedback loop 

first-order-servo treinsfer function 

integrating- servo transfer fxmction, Si(D) = i S(D) 

primary airplane transfer function 

numerator of G.(d) 

denominator of G(D) 

secondary airplane transfer function 

exponent of integrator in con^jensating network 

degrees of freedom of airplane 

externally applied disturbance function 

output of feedback-loqp compensating network 

servo time constant, seo 

compensating-network gain constant 

arbitrary polynomial defining airplane response, used in 
feedback- loop compensating network 

velocity-command gain constant, sec“^ 

acceleration-command gain constant, sec"^ 

wing span, ft 

wing area, sq ft 

relative density factor, -SL 

pSb 

mass of airplane, W/g, slugs 
weight, lb 

acceleration due to gravity, ft/sec^ 

nondimens ional radius of gyration in roll about longitu- 
dinal stability axis, 
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nondimens ional radius of gyration in yaw about vertical 
stability axis, ^K^q^cos^tj + KxQ2sin2T^ 

nondimens ional product of inertia, ” ^Z^^^sin tj cos tj 

nondimensional radius of gyration about principal X body 
axis of airplane 

nondimensional radius of gyration about principal Z body 
axis of airplane 

roll angle of airplane 

sideslip angle of airplane 

yaw angle of airplane 

angle of attack 

inclination of principal longitudinal axis with respect 
to flight path 

flight-path inclination from horizontal, 0 - a 
aileron deflection angle 
rolling angular velocity, radians/sec 
yawing angular velocity, radians/sec 
altitude, ft 

mass density of air, sl\jgs/cu ft 
true airspeed, ft/sec 
Mach number 

rolling-moment coefficient, moment 

qSb 

yawing-moment coefficient, 

qSb 

lateral-force coefficient, lateral force 

" qS 

trim lift coefficient 
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Subscripts: 

i 

0 

1 


input 

output 

limiting value of a variable 


A dot over a symbol Indicates differentiation with respect to time 
for exan^le^ fb = 


All angles are measured in reidians unless otherwise noted. 
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BASIC CONSIDERATIONS 


The automatic control system studied in this investigation incor- 
porates a compensating network, a device that eliminates airplane dynamics 
from the system response. The compensating network is essentially a 
computer that solves a set of eqiiations that are the inverse of the 
equations of motion of the airplane. The input to the network is either 
an gmgnii».r displacement, velocity, or acceleration and the output is the 
control-surface deflection required to give the input value of the 
variable . 

Figure 1 presents the block diagram of an automatic control system 
incorporating a forward-loop conpensating network. The dynamic response 
of this system is studied by a linear analysis and by simulation on an 
analog computer where nonlinearities and imperfect compensation were 
introduced into the problem. Inasmuch as the reg\ilatory response charac- 
teristics, the return of the airplane to a specified steady-state condi- 
tion when disturbed from that condition by an externally applied moment, 
may be of some interest, this type of control-system operation was studied 
in addition to the command-response dynamics. Appendix A presents the 
theoretical analysis of the forward-loop ccmpensating-network type of 
control system and the transfer function is developed. In addition, the 
feedback- loop compensating network is subjected to a brief analytical 
study in appendix A. This type of conpensating-network control system 
appears to eliminate some of the disadvantages of the forward-loop 
compensating network but mechsinization difficulties may bar it as a 
practical system. 


Application of the Ccmpensating Network to Roll Control Systems 

In order to evaluate a specific application of the conpensating 
network for the automatic control of airplanes, the principles set forth 
in appendix A were applied to three automatic roll control systems. 

These roll control systems differ in the number and types of feedback 
used and in this report the roll control systems are identified by the 
hipest order feedback used. Thus the basic control system (fig. 2(a)) 
is a displacement control or bank-angle -feedback system. The second 
system (fig. 2(b)), herein called the velocity control system, is 
derived from the first by the addition of a rolling-velocity feedback 
and a gain Ki. The third system (fig. 2(c)), called the acceleration 
control system, was obtained by adding a rolling-acceleration feedback 
and a gain K 2 to the velocity control system. 

In a conpensating network control-system feedbacks are not required, 
as in other types of control systems, to add damping and thereby improve 
the system response, since the cancellation of the airplane dynamics 
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eliminates the problems of poor airplane stability anddan^)lng. Thus, 
feedbacks are always made to the input side of the compensating network, 
since feedbacks to the servo for additional damping are unnecessary and 
if made would destroy the cancellation of the airplane dyneunics. The 
feedbacks shown in figure 2 are used to determine the airplane variable - 
bank angle, rolling velocity, or rolling acceleration - that is used in 
the compensating network to compute the desired control-siirface deflec- 
tion. This value of 6 is determined so that airplane dyneunics are 
eliminated from the system response as the airplane assumes the value 
of the variable that is fed to the compensating network. 

Since these roll control systems are considered as a part^ of a 
conqplete interceptor system, it is desirable to correlate the simulation 
used for the roll control systems with the automatic interceptor system. 

In an actual interceptor system the outputs of the airborne radar and 
director computer determine an azimuth error and an elevation error. 

These errors are used in the command ccmqjuter to define a roll command 
and a normal-acceleration command. The roll command depends upon the 
particular type of guidance system used and may be either a desired bank 
angle in space coordinates or a bank-angle error in interceptor coordi- 
nates. If an interceptor system is designed to command a bank angle in 
space coordinates, the bank-angle feedback is a necessary part of the 
roll control system. Therefore, in order to simulate this system, the 
reidar and computers need only be replaced by a bank-angle input to 
obtain a simple analog of the more complex system. However, if the 
radar-computer sy;stem operates in interceptor coordinates, the system is 
closed by a feedback loop through the radar. In this case the basic 
command to the roll control system is a rolling-velocity command that 
is proportional to the azimuth and elevation errors. Therefore, replacing 
the radar and con5)uter by a bank-angle input, a bemk-angle feedback, and 
a gain to convert the bank-angle error to a rolling-velocity command 
constitutes the. first approximation to the analog of the more complex 
system. Thus, the roll control systems shown in figure 2 are applicable 
to either one of the two guidance systems mentioned above. 


Simulator Setup 

Airplane equations of motion .- All airplane treuisfer functions, and 
consequently the transfer function of the compensating network, used in 
this investigation were derived from the linear equations of lateral 
motion with 7 = 0, which are as follows : 
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The general form of the airplane transfer function G(D) is 




G(D) = = aip2 + a2P + a^ 


^a(O) + b^D + be 


( 2 ) 


In equation (2) the coefficients a^ to a^ and b^^ to b^ are fvinc- 
tions of the mass and aerodynamic characteristics of the airplane and 
can be obtained by expanding the determinant of the airplane equations of 
motion, which is obtained from eqioations (l). 

Flight condition A, presented in table I, is the standard flight 
condition used in both the airframe and compensating-network equations 
when perfect compensation was desired. Incomplete compensation was 
introduced, aerodynamically, by substituting the singular or group varia- 
tions (flight condition A-1 or A-2) into the airplane equations. When 
it was desired to use two fli^t conditions in the problem, flight 
condition A remained in the conq)ensating network and flight condition B 
of table I was used in the airplane equations of motion. 

Compensating network .- The generalized transfer function of the 
conq)ensating network is 

F(D) = — - — ( 5 ) 

D«G(D) 

where G(d) is the transfer fxaiction ^/5a of the airplane (eq. (2)), 

K is the gain or an^lification through the network, 1/D is an inte- 
grator, and the exponent n prescribes the order of integration. In 
general, the maximum value of n that can be used is one greater than 
the highest order of the derivative of the bank angle that is fed back 
to form the input to the compensating network (see fig. 2). 
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The servo.- Thp servo vised throvigh the investigation is a first- 
order~time^lag servo, defined "by eqviation (Al) . The servo time constant 
is one of the variables of the investigation. 

Limiting was applied to the servo to simulate control-surface rate 
and displacement restrictions of a physical system. Two types of limiters 
were used, the winding and nonwinding types of limiters. When 6, the 
control-surface angular displacement, does not reach its limit, both 
types of limiters operate in an identical manner, even though &, the 
control-surface angular rate, is limited. However, when 8 reaches its 
limit, important differences occur in the operation of these limiters. 

In the winding- type limiter, when 6 > 5j the following condition 

J 8 dT < 8j (M 

must be satisfied before 8 can move off the stop. In the nonwinding 
type of limiter for S ^ &2 



8 = 0 


which implies that for 8 < 0, 8 moves off the stop immediately. It 

should be noted that the above discussion applies to positive limits. 

An equivalent representation can be written for the negative limits. 

These two types of limiters were included in the study because it was 
felt that winding limiters more closeiy approximate the operating charac- 
teristics of the proportional servo and its stroking motor when they are 
operating in a saturated condition, whereas the nonwinding limiter repre- 
sents a perfect proportional servo and stroking motor. 

nonwinding-type limiter as used in this problem was set up so 
that 8 goes to zero when 8 is limited. This retvirn to zero by 8 
was not carried throu^ to the recorders; consequently, 8 always shows 
a value on the record whether or not 8 is limited. Therefore, in order 
to determine . when 8 is limited from these records for the nonwinding 
limiter, the behavior of 8 must be taken into account as well as the 
behavior of 8. 
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SCOPE 


The investigation included both a linear analysis and a REAC study 
of the three types of roll control systems. The main purpose of the 
linear analysis is to show the effect of the system gains and servo time 
constant on the roll response of the airplane to a command input. The 
REAC study was conducted to investigate the effect of incon^jlete conden- 
sation and nonlinearities. 

Incomplete compensation results when the transfer function of the 
condensatlng network is not the exact inverse of the airplane transfer 
function. In the REAC simulation, incomplete compensation was introduced 
by changing the airplane flight condition from the flight condition for 
which the compensating network was designed and also by varying the 
stability derivatives singly and in combination in the airplane transfer 
function as shown in fli^t conditions A-1 and A-2 of table I. The 
latter simulates the problem which may arise in practice, where the 
condensating-network design is based on estimated stability derivatives, 
and the resulting network transfer function is not the exact inverse of 
the actual airplane transfer f\anction. The nonlinearities incorporated 
in the system are limits imposed on the maximum values of control-siarface 
rate and displacement. 


RESULTS AND DISCIKSION 
General 


The transfer functions used in the linear analysis of these control 
systems were derived as outlined in appendix A, assuming perfect compen- 
sation and linear operating conditions. The transfer function for each 
of these systems is as follows: 

For the displacement control system 


For the velocity 


^ K 

^ rp(l + tD) + K 

control system 


KKi 

h D^(l + tD) + KD + KKi 


( 6 ) 


(T) 
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For the acceleration control system 

^ ^ ^ 1 ^ ( 8 ) 

h rP(l + tD) + Kd 2 + KK2D + K1K2K 


Routh's discriminant was applied to the characteristic equation of 
each of these systems. Unstable conditions were indicated as follows; 

For the displacement control system 

n > 1 (9) 


For the velocity control system 

n > 2 

Kl< I 

For the acceleration control system 

n > 5 

K< Kit - 1 
Ki 

K< ^ 


,(n = 2) 


( 10 ) 


(n= 2) 
(n = 5) 


( 11 ) 


1 - tK2 

In addition, ■ system stability requires that all gains be positive. 


The rise time (the time required to reach 90 percent of the steady- 
state value) and the response time (the time required for the motions to 
reach and remain within 5 percent of the steady-state value) were used 
to analyze system responses and to make comparisons of different operating 
conditions. The linear analysis was used to determine the effects of 
varying the gains and servo time constant on the command response of the 
system. The variations were checked on the REAC and it was foimd that, 
when no limiting was present, there was good agreement between the linear 
and REAC results. However, when n, the exponent of the comgpensating 
network, is equal to 1 and a step input is used> initial conditions arise 
in the problem that destroy compensation. This situation, with REAC 
records, is discussed in the section on the displacement control system, 
and an analysis of this condition for the same control system is presented 
in appendix B. When limiting of the control^surface rate and displacement 
was added, the problem became nonlineeu:. The REAC was used to investi- 
gate this phase of the problem. Therefore, only REAC records showing the 
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effect of varying the limits on the control-surface rate and displacement 
are presented in this report. 

Equation (a 6) partially predicts the limiting condition introduced 
by increasing the gains because it shows that the forward-loop gain 
appears as a multiplicative factor in the an^jlitudes of 6 and 6. 

In general, the comments relative to the gains apply in the s«ime 
way when t, the servo time constant, is varied. However, REAC resvilts 
indicate that the value of t that is optimum in the linear analysis 
may not be the best when the control-s\irface rate and displacement are 
limited. For the displacement control system, where limiting is present, 
the oscillations introduced by t combine with the oscillations caused 
by the rate limiting to produce a system response that approaches the 
linear result. This effect is shown in figiire 3* 

When the roll response of the airplane is compensated, the ^ and 
'if motions are uncompensated. This occurs, as indicated by equation (a 6), 
becaxise airplane characteristics appear as a factor of the characteristic 
equation of the p and if transfer fimctions. 

Incomplete compensation was introduced by changing the stability 
derivatives individually or in groups or by changing the flight condition 
as indicated in table I. The resulting system response was in general 
unsatisfactory. However, in the case of the n = 5 acceleration control 
system, a satisfactory system response was obtained for the mixed-flight- 
condition type of incomplete conqpensation. More detailed resvilts of this 
study £ire presented in the sections on the displacement and acceleration 
control systems. The regulatory response for the basic control system 
and the two variations studied were uncompensated, as indicated by equa- 
tion (a 8) . Although the response was vmcompensated, it showed stable 
characteristics. For the magnitude of the disturbances considered 
(Cj = Cji = 0.1) limiting of the control-surface rate and displacement 
was troublesome only at hi^ -values of the forward-loop gain or servo 
time constant, or both. Since the regulatory response was uncon 5 )ensated, 
this mode of control-system operation was not extensively investigated 
and no results are presented. 

The linear analysis and the BEAC investigation of the effects of 
limiting are presented in the following sections for each of the control- 
system variations considered. When applicable, results on inccmqplete 
compensation are also presented. 


The Displacement Control System 

Linear analysis .- Equation (6) subject to condition (9) was used to 
determine the effect of -varying K and t on the response of the linear 
system. A step fvmction, = 60°, was used as the command input. 
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Figure 4 shows that, as K is increased from 5 to 12, both the rise 
time Sind the response decrease. As t is increased from 0.01 to 
0.5 second (fig. 5) > an oscillatory mode is introduced into the system 
response. An examination of the characteristic eqioation of the system 

indicates that these oscillations occur for t > -i. and that the fre- 

4K 

quency of the oscillations is proportional to — Jl - 4Kt. 

2t * 

REAP study .- In order to evaluate the effects of limiting and 
incomplete condensation, the control system was set up on an REAP type 
of analog condiuter. To check the REAP setup, a rim was made under 
linear conditions. A step input of - 60° was applied to the system 
and, inasmuch as no limiting occurs, the 0 motion shown in figjure 6(a) 
should be directly comparable with the linear curve in figure 6(b) . A 
conqparison of these two curves shows that the frequencies and rise times 
of the two motions are radically different. An examination of the equa- 
tions of motion of this control system indicated that initial conditions 
arising in the compensating network from the use of the step input were 
causing incond>lete conqjensation. The analysis presented in appendix B 
shows that, in order for complete compensation to take place, ^±(o) = 0 
The command-input network was changed so (ii(t) = - e-^'t) and 0c 

was set equal to 60°. Figure 7 shows the time histories of the recorded 
variables for the exponential input. A comparison of the 0 motion 
with that shown in figure 6(b) shows that the linear and REAP results 
are now in approximate agreement. 

As indicated by equation (a 6) , the motions y, 3, 5, and 6 due 

to 0i, shown in figure J, will be imcon^iensated and contain the oscil- 
latory characteristics determined from (tE)2 + D + k) (ail)2 + a 2 D + a^) 
which is the characteristic equation of the system. Thus, in addition 
to oscillatory modes that come from the control system, oscillatory 
characteristics determined by aiD^ + a 2 D + aj, the numerator of G(D) , 
will also be present in the f and 3 motions. In terms of the mass 
and aerodynamic characteristics of the airplane, the coefficients ax, 
a 2 , and aj are 
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Por the airplsme and control system studied herein^ the aiirplane induced 
oscillations are predominant. The period of the airplane oscillation, 
for the fli^t condition considered, is I .65 seconds, which is very close 
to the period of the Dutch roll oscillation. 

The effect of limiting 6 and 5: The effect of limiting 6 and 

8 was studied by using the nonwinding type of limiter . Conqparing 
figure 7 where 6 is unlimited with figures 8(a) and 8(b) where 
is set at 100 deg/ sec and 50 deg/ sec, respectively, shows that the effect 
of limiting 8 is to introduce a low initial peak and an oscillation in 
the 0 motion. Further reductions in 8j cause more limiting to occur, 
which in turn causes the severity of the limiting-induced oscillation to 
increase . 

When 8 is limited and 8j held constant, the effects are the same 
as noted for limiting of 8. As 8^ is reduced from 20° to 10° (figs. 9(a) 
and 9(b)), the amplitude of the oscillations increases. Setting 8j 
equal to 5° (see fig. 9(c)) causes the entire pattern of the oscillations 
to change. The 6 motion now approximates a square wave. The high- 
frequency oscillation in the 0 and 'f> responses is the third harmonic 
of -t;he basic frequency of the 5 square wave, the second term of the 
Fourier series for a square wave , and disappears as soon as the 8 square 
wave decays. A reduction in 8j to 40 deg/ sec (fig. 9(^)) produced no 
major change in the motion observed in figure 9(c) > where = 5° and 
8j = 100 deg/sec. Thus, it appears that, for small values of 8j, the 
effect of reducing is not critical. 

Incomplete compensation: The aerodynamic stability derivatives of 

the aiiplane were varied singly and as a group as shown in table I. Only 
two of the derivatives, Cnp and Cj^, produced a noticeable effect on 

the system response, the effect of Cnp being more pronounced than the 
effect of Cxp- When Cnp increased to 0.52 (fig. 10(a)), a very 

li^tly damped himting oscillation was introduced in the j6 and 0 
responses. The characteristics of this oscillation and its effect on 
the response time are such that the system response is considered 
unsatisfactory. 

Decreasing Cnp to 0.25 or increasing C^p to -O.O 8 had the same 

effect on the system response. Figure 10(b) shows an example of this 
effect on the system response for Cjp = -O.O 8 . This oscillation is 

less persistent than the one introduced by increasing Cnp (^*^ 8 * 10 (a))- 
However, comparing the system responses for CXp = - 0.08 with a similar 
case for Cxp = -O.IO 6 (fig. 9(a)) shows that the oscillations in the 
0 and 0 motions are more li^tly damped for the incomplete-rconpensation 
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case (fig. 10(b)). Decreasing to -0.15 produces a negligible 

effect on the system response. 

When the flight condition of the airplane was changed from flight 
condition A to flight condition B (table l) , the oscillation shown in 
figure 10(c) occurred. The motion is completely unsatisfactory with 
6 and 5 being limited most of the time. 


The Velocity Command System 

linear analysis of the n = 1 system .- A linear analysis of the 
n = 1 system was made by using equation (7) subject to condition (lO) 
to determine the effects of varying the gains K and and the servo 

time constant t on the system response. A step input of ^i(t) = 6o° 
was used in this analysis. 

The effect of varying K and K]_ is shown in fig\ires U(a) and 
11 (b) for T = 0.01 and 0.3 second, respectively. The effect of 
increasing K is to decrease the rise and response times for small 
values of Ki at both values of t. This effect is much less marked 
at high values of Ki. For combinations of K, K^, and t which 
result in oscillatory motions, that is, when 4 KK]^t > (K + l)^, 
increasing K adds dancing to the system. 

The most pronounced effect of increasing Ki is to reduce the rise 
and response times. As can be seen from figxire 2(b), the gain Ki 
appears as a multiplicative factor on the bank-angle error, and thus 
the rolling-velocity command which controls the speed of response is 
directly proportional to K^. 

In order to improve the system response, Ki shoiild be increased 
to control the rise and response times while K is increased to add 
damping to the system. The servo time constant t should be kept as 
small as possible. 

REAC stTJdy of the n = 1 system .- The n = 1 system was studied 
on the REAC by using the nonwinding type of limiter and the exponential 
input . 

With 5j set at 20°, was decreased from 100 deg/sec 

(fig. 12(a)) to ho deg/ sec (fig. 12(b)). This reduction in Bj causes 
the system to become more oscillatory, and the rise and response times 
£ire increased. A comparison of figures 12(a) and 12fc) shows the effect 
of reducing Bj to 5° while remains at 100 deg/sec. For this 

condition the B motion closely approximates a square wave, and again 
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the third harmonic of the basic frequency of this wave appears in the 
^ and ^ motions. The rise time of the system is increased while the 
response time is almost xmchanged. Reducing from 20° to 5° when 

is deg/ sec (figs. 12(h) and 12(d)) produces results similar to 

those previously obtained. A comparison of figures 12(c) and 12(d) 
indicates that a reduction in when 6^ is 5° does not appreciably 

affect the system response. The square-wave phenomenon and the insen- 
sitivity to reductions in for = 5° were also noted in the dis- 

placement control system. 

Linear analysis of the n = 2 system .- Equation (7) subject to 
condition (10) was used to make a linear analysis of the velocity control 
system for n = 2. A step input of /6i(t) = 60° was used as the forcing 
function. As in the linear euialysis of the n = 1 system, K, and 

T were varied to determine their effect on the system response. Fig- 
ure 15(a) shows the effect of increasing t from 0.01 to 0.5 second on 
the system response. As t is increased the system response becomes 
oscillatory, a condition that would be expected since the condition for 

neutral oscillatory stability is ^ = K]^. Thus for constant the 

severity of the oscillation increases directly with t. 

As Ki is increased from I.5 to I6 (fig. 15(b)), the system 
response becomes oscillatory. However, in spite of the oscillatory 
mode, both the rise time and response time show large decreases. The 
oscillations may be reduced or eliminated by Increasing K. With 
Ki = 16 increasing K to 64 gives a deadbeat response for the system 
and further decreases the rise and response times. The gains K and 
perform the same functions for this system as they did for the 
n = 1 velocity control system, controlling the speed of response 

and K the damping. However, in the n = 2 system, in order to mini- 
mize the oscillatory characteristics in the response, when t is 
constant, K must be greater than Kj. 

REAP investigation of the n = 2 system .- The n =» 2 velocity 
control system was studied on the REAP with both types of limiters; a 
step input of }6i(t) = 60° was used as a forcing function. 

Two types of limiting were used for the n = 2 velocity control 
systeia, the nonwinding and winding types of limiters. These limiters 
have been discussed in a previous section of this paper. 

When 5 is made unlimited by setting t = 0, the two types of 
limiters are equivalent. This condition is used to show the effect of 
limiting B. Figures l4(a) and l4(b) show that, as 62 is decreased 
from 20° to 10°, the smooth response of figure l4(a) is modified by an 
oscillatory mode. A decrease in B2 to 5° (fig. 14 (c)), causes this 
oscillation to tecome more severe. 
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The effect of reductions in 6^ is shown for = 20° in fig- 
ures 15(a) and 15 (t) and for = 5° in 15(c) and 15(d) for the non- 
winding l^e of limiter. For 8j = 20° a reduction in Sj from 
100 deg/ sec to 40 deg/ sec causes oscillations to develop in the (6 and 
^ motions. When a similar reduction is made for = 5°# 00 appre- 
ciable change occurs in the motion. However, the square-wave condition 
at low values of 5^ noted for the preceding systems occurs in this 
case. The results for the winding type of limiter are shown in fig- 
ures 16(a) and (h) for = 20° and in figures l6(c) and l6(d) for 
62 = 10°. When bi is reduced from 100 deg/sec to l)-0 deg/sec for 
6j = 20°, the oscillations occurring in the ^6 and p motions become 
more severe. . When 5j is 10°, a reduction of 6 2 from 100 deg/sec to 
60 deg/ sec causes the system to become unstable. No stable cases for 
62 = 5° were foimd for the winding-type limiter. Until the system 
became unstable, the rise times varied very little, but, becaiise of the 
oscillations, the response ti™e was gradually increasing. 


The Acceleration Control System 

Linear analysis of the n = 5 system . - A linear analysis of n = 3 
acceleration control system was made by using equation (8) subject to 
condition (ll) . The system constants that were varied dirring this inves- 
tigation are K, Ki, K2, and t. 

Figure 17 shows the effect of varying K on the ^ response. As 
K is increased from 4 to 20, the oscillations and overshoot are elim- 
inated. This increase in K is accon5>anied by an increase in rise time 
and a decrease in response time. The oscillations that occur for K = 4 
would increase with decreasing K since for the values of Kg, and 

T used the system is approaching the condition for neutral oscillatory 
stability, which occurs at K = 1.55* 

As Ki is increased from 0.75 to 2.0 (fig. I8) , both the rise time 
and response time decrease. At = 2.0, an overshoot and sli^t oscil- 
lation occxrr in the response. Further increases in would cause this 

oscillation to become more severe. 

As Kg is increased from 5 to 6 (fig. 19) ^ the response gradually 
flattens vmtil at K2 = 6 a visible change in slope occurs at t = 0.8 
second. This change in slope is probably caused by a change in the danping 
introduced by the change in Kg and the limits placed on the velocity 
and acceleration commands because of the way in which they are determined. 
With Kg = 6, increasing Kx to 2 improves the system response suid 
decreases the rise and response times. The increase in Kx has decreased 
the system damping eind liberalized the limits on the velocity and acceler- 
ation commands . 
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As T is increased from 0.1 to 0.2 second (see fig. 20) , the 
response is practically the same up to some point between t = 0.1 to 
0.2 second. At this point an oscillation starts to develop, which is 
quite evident at t = 0.2 second (see dashed curve in fig. 20). As 
T is further increased, this oscillation will become more severe until, 
at T = 0.5 second, neutral stability occurs (see eq. (8) and conditions 
for \anstable operation) . 

A study of figures Y{, l8, and 19 indicates that K controls the 
damping, while K2 controls the rapidity of the initial response and 
affects the damping. The gain exerts an influence on the rise and 

response times that is inversely proportional to magnitude. Thus, 
increases in Ki and K2 will in5>rove the speed of response and 
Increasing K together with the increase in K2 keeps the system 
response nonoscillatory and should give a better overall system response. 
Accordingly, Ki was increased to 5*25 > K2 to 10, and K to 20 while 
T was held constant. This response, shown in figure 21, is deadbeat 
with much smaller rise and response times than have previously been 
obtained for the n = 5 acceleration control system. 

REAP investigation of the n = 5 system .- A step input of 
0i(t) = 60° was used as the forcing function in the REAC study of the 
n = 5 acceleration control system. A winding-type limiter was used 
for this system. Since the winding-type limiter gives a more conserva- 
tive result for a given system than the nonwinding type of limiter, the 
trends for the n = 5 system as studied should also indicate the trends 
to be expected if the nonwinding limiter had been used. The major 
difference between the effect of the two types of limiters is that 
regions of unstable operation indicated for the system with the winding- 
type limiter might not occur vmder similar conditions for a system 
equipped with the nonwinding limiter. 

The effect of limiting 5 and 6: In order to show the effect of 

reducing on the system response, the limit on the control-surface 

rate was removed by setting x = 0. As is reduced from 20° to 10° 

(figs. 22(a) and 22(b)), no apparent change occ^s in the ft motion 
and a very slight oscillation is noted in the p and f) responses. 

When 6 2 was reduced to the system became unstable. 

With hi set at 100 deg/ sec and = 20° (fig- 25(a)), a non- 
oscillatory response that is almost equivalent to the bj = <» response 
(fig. 22(a)) was obtained. Reducing to 40 deg/sec (fig. 25(b)) 

introduces an extremely small oscillation into the and motions. 

Available results indicate that, for = 10° and = 100 deg/ sec 
and 40 deg/sec, the system response does not change appreciably from 
that shown in figure 25. When was reduced to 5°, unstable condi- 

tions were encountered for all values of 6^. 
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Incomplete compensation: When the stability derivatives were vsoried 

individtially or in a group, as indicated in table I, the n = 5 acceler- 
ation control system response was modified in the same way as the response 
of the displacement control system. 

Changing the airframe flight condition from A to B (table l) 
produced unstable resxilts for low values of K, thus following previous 
trends. However, when K was increased to l4 (fig. 24(a)), the insta- 
bility was reduced to a hunting oscillation. Increasing K to 24 
(fig. 24(b)) practically eliminated these oscillations in the 0 response, 
and the oscillation in the other recorded variables is damped. 

Remarks on other control systems using acceleration feedback . - 
Acceleration control systems described by equation (8) with n = 1 or 2, 
were superficially investigated by a linear analysis. The general trend 
of these reduced-order acceleration control systems was to reqvilre a 
different distribution of the forward-loop gain between K^, K 2 , and K 

than for the n = 5 system for the most satisfactory response. In 
addition, there are indications that higher values of the forward-loop 
gain can be used before limiting oscillations become severe enough to 
affect the system response adversely. 

An acceleration control system with the ^ feedback eliminated and 
K 2 = 1 was investigated by linear methods. Routh's criterion indicated 
unstable conditions for n > 1. For n = 1 the system was stable for 
very small values of K and the allowable reinge of variation of K was 
small. The response of the system was very slow unless extremely large 
values of Ki were used. These large values of Ki commanded airplane 
accelerations that are considered as too high to be practical. Because 
of the high accelerations and the narrow band of K for stable operation, 
this system was not investigated on the REAC. 


CONCLUDING REMARKS 


The concept of the compensating network, a computing device that 
eliminates the dynamics of the controlled element from the system response, 
has been applied to automatic control systems. The general characteristics 
of this type of control system have been determined and the results applied 
to the analysis of three related automatic roll control systems for air- 
planes. These roll control systems differ in the number and type of feed- 
backs that are used to svqpply information to the compensating network for 
use in computing the control orders tha't produce cancellation of the air- 
plane dynamics. 
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Within the limitations imposed by control-surface rate and dis- 
placement limiting and by imperfect coD 5 )ensation, the compensating net- 
work can give adequate command roll control that is free of airplane 
dynamic characteristics. Limiting of the control-surface rate and dis- 
placement introduces oscillatory modes into the system response that 
become more severe as the limiting time increases. Because of the way 
in which the limiting affects the system, restrictions were placed on 
the gains, the order of integration in the compensating network, and the 
servo time constant. The most critical effect caused by imperfect com- 
pensation occiurs when a change in airplane fli^t condition, not accounted 
for by the compensating network, introduces an unstable response. Inaccu- 
racies in the airplane mass and aerodynamic parameters used to design the 
network are not critical except in the case of Cnp and 0 %^ (partial 

derivatives of yawing-moment and rolling-moment coefficients with respect 
to sideslip angle) when hunting oscillations are introduced into the 
response. In addition to the difficulties introduced by limiting and 
incomplete compensation, the control of the yaw and sideslip motions are 
characteristically uncompensated in the command mode of operation. 

The effects of limiting and imperfect con 5 >ensation on the regulatory 
response are the same as for the command response. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 12, 1955- 
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APPENDIX A 

DEVELOPMENT OF TRANSFER FUNCTIONS FOR COMPENSATING-NETWORK 

CONTROL SYSTEMS 


In this appendix the transfer functions of automatic control sys- 
tems using compensating networks are developed. The transfer functions 
are derived for linear operating conditions and perfect compensation. 
The forward- and feedback-loop coD 5 )ensatlng network are discussed in 
the order mentioned. 


The Forward-Loop Compensating Network 

The block diagram of an airplane equipped with an automatic con- 
trol system is shown in figure 1 where F(D), S(d), and G(d) are 

the transfer functions of the compensating network, servo, and airplane, 
respectively. 

An examination of the response characteristics of typical high- 
performance servos Indicated that over the range of airplane frequencies 
a first-order time-lag servo was a good first approximation of a phys- 
ical servo. Accordingly, the servo transfer function was taken as 


S(D) 


1 

1 + tD 


(Al) 


where t is the servo time constant. The closed-loop transfer function 
of the system is 


^ ^ F(D) S(D) G(D) 
Xi 1 + f(D) S(D) G(D) 


(A2) 


It is desired to determine the transfer function of the compensating 
network F(D) so that the airframe dynamics are canceled and that as 
a command system the closed-loop response has zero steady-state error. 
Basic servomechanism theory requires that, in a closed- loop system, an 
integration take place to satisfy the zero-steady-state-error condition 
(see ref. 5)» Therefore, a logical choice for the conipensating-network 
transfer function is 


K 

dg(d) 


F(D) = 


(A5) 
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Svibstituting for F(d) and S(D) in equation (A2) results in the 
closed-loop transfer fxmction 



K 

tD^ + D + K 


(Al^) 


and 


lim ^ = 1 (A5) 

D— >0 % 


Th\is, the requirements that the airplane transfer is canceled and the 
steady-state error is zero are both satisfied. It should be noted that 
in corapensating-network control systems the integration must always be 
introduced explicitly because implicit integration in the airplane trans- 
fer function is eliminated by the cancellation process. 

In the foregoing analysis the transfer function F(d) was deter- 
mined to give complete compensation for a particvilar degree of freedom 
of the airplane X. If there is another degree of freedom t), which is 
related to X by the transfer function Gq(D), (see fig. l), ’lo/^i 

given by 


Hq _ 

Xi ■ Xi Xq 


Gi(D) = 


K 


Gi(D) 


+ D + K 


(a6) 


If Gq^(d) is the ratio of two polynomials, the output motion t)q 
will contain modes determined by the denominator of Gi(D) since it 

appears as a factor of the characteristic equation of the system, and 
thvis the response is uncompensated. 

Bebn-vior of the compensating-network control system as a regulator. - 
Tn addition to providing compensated control in response to command inputs, 
the control system is sometimes called upon to act as a regulator. If a 
distvirbance M(D) is applied to the airframe, as shown in figure 1, and 
the input Xj = 0, the closed-loop transfer function is 

G(D) 


M 


1 + f(d) s(d) g(d) 


(A7) 
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Substituting S(d) = and F(d) 

1 + tD 


K 

DG(D) 


into equation (A7) yields 


^2 = . D(l . :+ tD) 

^ td2 + D + K 


(A8) 


If G(D) is the ratio of two polynomials in D and the denominator 
contains a constant, the steady-state condition for the regular response, 

Y 

lim ^ = 0, is satisfied. However, G(D) that appears in equa- 
D-^O ^ 

tion (a 8) is a transfer function of the airplane. Since the denominator 
of G(D) is a factor of the characteristic equation of the system, the 
response will contain airplane characteristics. Thus, the forward-loop 
compensating network does not eliminate airplane dynamics from the regu- 
lator response. 

When the control system is responding to external distiurbances 
(that is, when moments are applied to the airplane) the motion tJq is 

related to the disturbance M by the following transfer function: 


^O _ ^ Ho - ^ P (Vi) 
M M M ^ 


G(D) Gi(D) 

+ D + K 


(A9) 


Thus, motions determined in part by the denominator of G(D) will 
appear in the t 1 q(d) motion. Thus, in most cases the response 

will contain airplane modes and therefore the response is uncompensated. 


The Feedback-Loop Compensating Network 

Figure 25 is a block diagram of a simple displacement type of con- 
trol system incorporating a feedback-loop compensating network. The 
transfer function of this control system for the command mode of opera- 
tion, 0, M = 0, is 


^2 = Si(D) g(D) 

% 1 + Fi(D) Si(D) G(D) 


(AlO) 


and for the regulatory mode of operation, Xj^ = 0, M ^ 0, the transfer 
fvinction is 
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^ (All) 

” 1 + Fi(D) Si(D) G(D) 

where S^Cd), G(d), and F 3 ^(d) are the transfer functions of the servo, 
airframe, and compensating network, respectively. 

For the feedback-loop compensating-network control system, the 
servo transfer function assumed for this study is 

Si(D) = i S(D) = (A12) 

D D(1 + tD) 


and if the transfer function of the cornpensating network is assxamed to be 

, R(D) s,(D) - mw 
KSi(D) N(D) 

the following system transfer functions result; 


^ _ kh(d) 

Xi R(D) 

r ' ^ ^ ^ gi(d) 

Xj_ Xq 


for the command responses and 


3^ _ KD(1 + tD) N(d) 
M r(D) 


(Aik) 


(A15) 


(Al6) 


To _ ^ To ^ KD(1 + tD) N(D) / V 
M M Xq r(D) ^ 


(AIT) 


for the regulatory responses. 

In these transfer functions R(d) is an arbitrary polynomial in 
D that defines the airplane response. If R(d) is taken so that 
(l + tD) appears as a factor, the transfer functions for the regulatory 
response will become simpler. 
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Perfect compensation is present in contrast to the lack of com- 
pensation for some phases of operation found in the forward-loop com- 
pensating network obtained for all these responses. For equations (a 1^) 
and (A16) there are no airplane characteristics in the characteristic 
equation of the system. In equations (A15) and (AI 7 ) the denominator 
of Gi(D), which is N(D), is canceled by n(D) in the numerator. In 
addition, equations (A11<-) to (Al?) meet the required steady-state condi- 
tions for the command and regulatory responses. 

It would appear that the feedback-loop compensating network^is-an" 
ideal compensating system. However, there are two practical difficulties 
that arise with respect to this system. First, it appears that R(D) 
must be a seventh-order polynomial in order to prevent derivatives from 
appearing in the system transfer function. Second, in the compensating 
network the numerator has a hi^er order than the denominator, which 
introduces derivatives into the system and the order of the derivative 
is such that it might make it very difficult to mechanize the compensating 
network. The use of a seventh-order polynomial for R(D) wovild aggravate 
this condition. 

Thus, before the feedback-loop compensating network can be evalu- 
ated, more basic research is required to determine whether the above- 
noted difficulties are inherent in the system. 
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APPENDIX B 

THE INITIAL CONDITIONS IN THE COMPENSATING-NETWORK 
CONTROL SrSTEM 


In the REAC study of the displacement control system, it was found 
that when a step input was used the response of the control, in the 
absence of control -surface rate and displacement limiting, was incon- 
sistent with the results of the linear analysis for a similar input. 
Since servomechanism theory prescribes that, in order for a transfer 
function, as used in this paper, to exist, all initial conditions in 
the system must be zero (see ref. p. 84), the characteristics of the 
step input and the equations of motion of the system were examined to 
determine if the assumption of zero initial conditions through the sys- 
tem had been violated on the REAC. 

The step input was examined (see ref. 6) and it was foxmd that for 
t = 0 the step input has an average mean value of l/2. This means that 
for the 6o° step used the value of at t = 0 is 50°. 

Figure 26 is a block diagram of the displacement control system as 
it was set up on the REAC. The equations of motion of this system are 
as follows: 


1 

II 

(Bl) 

£i = Yi(^^ - 5Zlo) 

(B2) 

Vi = K(jf5i - - Xo 

(B5) 

Vo = iig(D)[K(9(i - ^o) - 

(b4) 

(d - ki)Xo - (kgD^ + k^D) 1 - k^^ = k^Vo 

(B5a) 

-(kgD + k^)Xo + (d 2 - k8D)| - k^^ = 0 

(B5b) 
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-^loXo " 1^11 ^ Xq - - 0 


(B5c) 


Vq = (1 + TD)&a 

- (^2^^ + kjD^t - = k^B^ 

-(k^D^ + kyD)^ + (l>2 - kgD)if - k^p = 0 
■ (^ICp + ^u) ^ pD ^ + - k]^ 3 = 0 


(b6) 


(B7) 


where k]^ to k]Lj are fvinctions of the airplane mass and aerodynamic 
parameters, e, , and Vq are defined in figure 26 , and Xq, |, 

and 5 correspond to \jf, and p, respectively. The transfer func- 

tion of the high-gain amplifier is g(D) and li is the gain of this 
amplifier. 

After assuming the transfer function of the hi^-gain amplifier g(D) 
to be unity and taking account of the initial condition on the deter- 

minant of the equations of motion was expanded for ^q(d). This expan- 
sion gave 


9*0 = 


-fl(D)|^lK[J f2(D) + tifi(D)] - n2Kfi(D)|;Zii 


f2(D) 


-(1 + tD)|| f2(D) + nfi(D^ - fi(D).nK[| f2(D) + nfi(D^ - 


-jifl(D). 

[-(k2D - k5)(D - 

ki4) - k4ki2 

k6.^ri(l» 

•Xo(0) 

f2(D)' 

-(1 + tD)[| 

f2(D) + ufi(D^ 

1 - fi(D)^ 

.k[ 

^ f2(D) + Mfi(D) 

- ^2KfJ^(D) . 


• (b8) 
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and 


lim 




tD^ + D + K 


+ 





(kgD - kj)(D - k^i^) - ki^ki2_ 

fl(O) 

■Xo(0) 


(tD2 + D + K)f2(D) 

(B9) 


In this equation fi(D) is the numerator of the transfer function 
g ntl f2(l>) is the denominator of that transfer function. 


The first term on the right-hand side, after being divided by 
is easily recognizable as the transfer function for 0o/^i used in the 
linear analysis. The second term is a fimction of the initial value 
of Xq. If Xo(0) = 0, this term is zero and complete compensation takes 

place. If Xo(o) ^ 0, this term is not zero and complete compensation 
does not occxir. By evaluating equations (B1) to. (B5) at t = 0, it can 
be shown that 


Xo(0) = K0i(O) 


(BIO) 


Thus, if a step input is used, complete compensation will not take place. 

However, if J^i(t) = ^^(l - e-^t), then 0^(0) = 0 and from equa- 
tion (bio) Xq(o) = 0; thus the second term on the right-hand side of 
equation (B9) becomes zero and complete compensation takes place. 

This analysis was extended to cover the velocity and acceleration 
control systems. For these systems it was found that when n = 1 the 
exponential input had to be used to get complete compensation. For ^ 
n > 1 the step input co\ild be used and complete compensation was / 
obtained. 
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TABIE I.- AIRPLAME FUGHT CONDITIONS AND 
CHARACTERISTICS CONSIDERED 


Parameter 

Fli^t condition 

A 

(standard) 

A-1 

(singular 

variations) 

A-2 

(group 

variations) 

B 

H, ft 

60,000 



60,000 

p, slugs/cu ft 

0.000224 



0.000224 

W, Ib 

26,547 



26,547 

b, ft 

55.81 



35.81 

S, sq ft 

4oi 



4oi 

V, ft/ sec 

1,942 



1,359 

e, deg 

6.95 


( 

0 

CL , deg 

5.58 



5.04 

q = 6 + a, deg 

, 10.53 



5.04 


256 




ro 

. 0.0151 



0.0124 


0.115 



0.116 

Kxz 

0.0188 



0.00916 

b/v 

0.01844 



0.02655 

M 

2.0 



1.4 

Clu 

0.044 



0.0655 

Cl 

0.157 



0.520 


-0.0775 



-0.117 


-0.205 

-0.15,-0.25 

-0.15 

-0.275 

% 

-0.106 

-0.08,-0.15 

-0.15 

-0.128 

CZr 

0.158 

0.12,0.2 

0.12 

0.189 

Cnp 

0.0275 

-0.01,0.05 

-0.01 

-0.014 

Cnp 

0.285 

0.25,0.52 

0.52 

0.545 

Cnr 

-0.600 

-0.5,-0.7 

-.500 

-0.690 

% 

-0.695 

-0.6,-0.8 

-0.600 

-0.785 
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Figure 1.- Block diagram of generalized automatic control system incor- 
porating a compensating network. 








(b) Velocity control system. 



(c) Acceleration control system. 


Figure 2.- Block diagrams of roll control systems incorporating a com- 
pensatory network. 














Figure 3»- Effect of t on the response of the displacement control system. 

5^ = 100 deg/sec; flight condition A; exponential input a = 
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eap‘ (f, 


flight condition A. 


t, secovids 


0 


t, secoyuis 

(a) REA.C result for step input = 60°. 

Figure 6.- Response of displacement control system. K = 5j t = 0.1 second 

flight condition A; no limiting. 



57 



6dp ' (j) 


Figure 6.- Concluded. 


- 2 . 
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e 7.- Response of displ 
= 8; flight condition A 

miting. 


“i^dea/sec cfjdej 




Figure 8.- Effect of limiting of 6 on the response of the displacement control system. 
T = 0.1 second; flight condition A; exponential input a = 4 . 




Figvire 9.- Effect of limiting of 6 on the response of the displacement control sysh 
T = 0.1 second; flight condition A; exponential input a = 6. 





Figure 9*- Concluded 
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10.- Effect of incomplete compensation on the response of the displacement control system. 
K = T = 0.1 second; exponential input a = 6; 6^ = 20^; 5^ = 100 deg/sec. 
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Figure 12.- Effect of limiting of 6 and & on the response of the n = 1 velocity control sys- 
tem. Ki = 4.0; K = 8.0; t = 0.1 second; flight condition A; exponential input a = 6. 
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(a) bi = 20°; 6j = 100 deg/sec. (b) 6j = 20°; 6j = ^t-O deg/sec. 

Figure 15. - Effect of limiting of 5 on the response of the n = 2 
velocity control system. = 1.5» K = 8.0; t = 0.01 second; fli^t 
condition A; step input (ft a = 6o°; nonwinding limiter. 
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Figure 17- - 0 response of the n = 5 acceleration control systems for 
various values of K. Ki = 1 . 5 ; K 2 = t = 0.01 second; flight 
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Figixre 19* - 0 response of the acceleration control system for various 
values of K 2 and Kx* K = 8.0j t = 0.01 second; fli^t condition 
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